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^ Abstract 

Assuming Majorana nature of neutrinos, we re-investigate, in the light of the 
recent measurement of the reactor mixing angle, the allowed ranges for the absolute 
values of the elements of the neutrino mass matrix in the basis where the charged- 
lepton mass matrix is diagonal. Apart from the derivation of upper and lower bounds 
on the values of the matrix elements, we also study their correlations. Moreover, 
we analyse the sensitivity of bounds and correlations to the global fit results of the 
neutrino oscillation parameters which are available in the literature. 
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1 Introduction 



Our knowledge of the neutrino oscillation parameters has enormously improved in the 
recent years. The experimental results of the Double Chooz, Daya Bay and RENO Col- 
laborations [1] have impressively confirmed the earlier hints [2] for a non-zero reactor 
mixing angle. Taking these novel results into account, the recent global fits [H H] estab- 
lish #i3 > at a confidence level of ~ 10<r. Also the values of the solar mixing angle and 
the mass-squared differences are now known with good accuracy. While for a long time 
sin 2 #23 ~ 1/2 has been a very good approximation for the best fit value, the recent global 
fits [31 H] hint towards a deviation from maximal atmospheric neutrino mixing. However, 
even at la it is not clear from the recent fits in which octant #23 lies. The global fit of [3] 
allows #23 = 45° at 2a and that of [I] allows maximal atmospheric mixing within the 3a 
range. The least known mixing parameter is the CP phase 6: it is unconstrained at 2a, 
but the best fit values of [I] hint towards 5 ~ 71. 

With all these improved results, it is worthwhile to perform an investigation of the 
allowed ranges for the elements of the neutrino mass matrix M. v , as was done by Merle and 
Rodejohann in their seminal paper [5]. In the context of textures of M. v , the correlations 
of the elements (M.v) a p (a, (3 = e,/i, r) of the neutrino mass matrix are of particular 
interest. Therefore, the main goal of our paper is the construction of plots correlating 
the absolute values of the elements of the neutrino mass matrix with each other. We will 
assume Majorana nature of neutrinos in this paper. 

The structure of our paper is as follows. In section [2] we investigate analytically upper 
and lower bounds on the absolute values of the elements of the neutrino mass matrix. In 
section [3] we will outline the numerical methods we will apply in our analysis, and the 
results will be presented in section HI 



2 The elements of the neutrino mass matrix 

Assuming neutrinos to be of Majorana nature, the neutrino mass matrix is a complex 
symmetric 3x3 matrix, which can be diagonalized as 

U T M U U = diag (m l , m 2 , m 3 ) , (1) 

where the m ; are the neutrino masses and U is a unitary matrix. In the following we will 
always assume to work in the basis in which the charged lepton mass matrix is given by 

Mi = diag(m e , m M , m T ), (2) 

which implies U = £/pmns- Consequently, U can be decomposed as 

U = D 1 VD 2 , (3) 

where D\ = diag(e l</3e , e lip ^, e tlfT ) is a diagonal phase matrix, V is the mixing matrix in 
the parameterization fixed in [8], and D2 = diag(e J<J1 , e lf72 , e lf73 ) is the matrix of Majorana 
phases. Without loss of generality, we assume 03 = 0. Inserting this into equation ([1]), 
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leads to 



(M v )ap = (^diag(m 1 , m 2 , m,)^) = e ~^^) ^ m fe e" 2 ^^^. (4) 



fc=i 



The absolute values of the elements of the neutrino mass matrix 
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(5) 



depend on nine real parameters, namely the three neutrino masses rrii, the three mixing 
angles 8 12 , 623 and #13, the Dirac CP phase 5 and the two Majorana phases o\ and Oi- 

From equation ([5]) we can deduce an upper bound on ((.M^)^! as follows. Rewriting 
the absolute value as a scalar product 



\(M 
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we can use Cauchy-Schwarz's inequality to find 

\{M v ) aP \<\A\\B\. 
Due to the unitarity of V, we have J2 k \ V ak \ 2 = 1, and thus 



\(A\B)\, (6) 



(7) 



\A\ < 4 /maxmfc, \B\ < /maxm^, 



which leads to the final result 



\{M v ) a R\ < maxm fc , 

k 



(9) 



i.e. the absolute value of an element of the neutrino mass matrix is smaller than the largest 
neutrino mass. 

We can also construct a lower bound on \{M. y ) a p\. Defining 



a k = m k \V ak \\Vp k \ 



(10) 



and taking into account that the Majorana phases are not constrained by experiment up 
to now, we find 
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where the p k are unconstrained phases. Now we have to distinguish two cases. If the three 
numbers a k are such that they can be conceived as the lengths of the sides of a triangle, 
then the right-hand side of equation fTTTT) can become zero and |CM„) Q £| as well. It is 
easy to show that it is possible to construct a triangle with side lengths a k , if and only i£] 



2 max a k — a k < 0. 



(12) 



1 A different but equivalent inequality has been derived in [B]. 
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Therefore, we end up with the 




\{M v ) a p\ > 2maxa fc - ^a fc . 



(13) 



k 



In this way, one gets rid of the Majorana phases and this inequality may be used to rule 
out single texture zeros in the neutrino mass matrix. 

3 Numerical analysis 

As already discussed in the previous section, the absolute values of the elements of M. v 
depend on nine variables, two of which — namely the two Majorana phases <j\ and o<i — are 
totally unconstrained. In [3] the then available experimental data were used to produce 
plots of the K-My)^! versus the smallest neutrino mass m . The goal of the present 
paper is to produce — using the results of the latest global fits of neutrino oscillation 
experiments — plots of ((-M^)^! versus \{M. v )a.'P'\- Since the knowledge of the neutrino 
oscillation parameters has improved considerably in the recent years, we also redo the 
numerical analysis of [5] and show the plots of K-Mj,)^! versus rrtQ. 

Let us now turn to our numerical strategy. Concerning the desired plots, a first attempt 
of creating scatterplots was, unfortunately, doomed to failure because even random point 
numbers as high as 10 9 were not sufficient to fathom enough of the allowed parameter 
space to achieve appealing plots. Therefore we follow a different strategy. From the 
scatterplots we can guess the shapes of the areas which would be filled in the limit of 
infinitely many points. In particular, we find that the allowed areas have no "holes," from 
where it becomes clear that it is sufficient to construct their boundaries. We do this in the 
following way. Consider a plot showing |(A^i/) a /3| or m on the x-axis and \{M.y)a'i3'\ on 
the y-axis. Then we start by pinning the quantity on the x-axis to some given value xq. 
Then we minimize and maximize y for fixed x = xq and obtain two points (xo,y m in) and 
(x , Vmax) of the boundary. Afterwards we do the same for fixed y = y which leads to the 
two points (x min ,y ) and (x max , y ). Repeating this procedure with a suitable number of 
different values for x and y$ finally yields the desired allowed area. 

As described above, we need an algorithm which allows to minimize and maximize 
real functions of nine variables^ For this purpose we choose the Nelder-Mead algorithm 
(downhill simplex method) [TJ, which is a direct search method for finding a local minimum 
of a given function. However, since all functions / we will consider are non-negative, by 
minimizing — / we can use the Nelder-Mead algorithm to find a local maximum of /. 
Since we are interested in the global minima (maxima) of /, single runs of the algorithm 
are not sufficient. Thus for every minimization (or maximization) we start with 200000 
different random start simplices and also perform perturbations of candidates for a good 
minimum (maximum). Since there are nine parameters in Ai u , the domain of the Nelder- 
Mead algorithm is, by construction, the full parameter space R 9 . In order to restrict the 

2 If the right-hand side of inequality ([T3f is negative, then the resulting bound is |(Afi/) a ^| > 0. 

3 Analysing the predictions of the "best fit" results of the global fits would reduce the number of 
parameters. However, instead of adapting our program to a smaller number of parameters, we let the 
parameters which shall assume their best fit values b vary in the interval (b — e\b\, b + e\b\), with e = 1CP 6 . 
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search to some domain D C E 9 , we add the characteristic function 




for p E D , . 

n J>(P : =<^ . P , n 14 
for p f. D 

to the function which has to be minimized (maximization of / then is equivalent to 
minimization of — /). The vector p collects the nine parameters of the | {M. v ) a p\. The 
pinning of x to Xq is also achieved by adding a characteristic function, namely 

Hz(x(p)) := K *>Tx(p) Gl 

A ^ 1Q 6 for x ^ g j> V J 

where 7 is an interval containing xq. The pinning of y to yo is done in the same way. 

We also want to evaluate the lower bound (fT3l) derived in section [2j For this purpose, 
we will use the Nelder-Mead algorithm on the function 

2 max a k - ^ a k + U D (p r ) (16) 

k 

of the seven real variables p' = (m , Am^, Am^, 6*12, #23, #13, 5). 

Next let us discuss the domain D of the parameters p. For the mass-squared differences 
and the sines squared of the mixing angles we use the best fit as well as the ncr-ranges 
(n = 1,2,3) provided by [H H]. The Majorana phases o\ and 02 are unconstrained and 
can thus vary between zero and 2tt. The situation for 5 is a bit ambiguous because the 
best fit values are very different in the two global fits [3] (version 3) and jl]; while for 
the normal ordering of the neutrino mass spectrum both fits favour a value of 5 which is 
roughly it, the best fit values in the case of an inverted spectrum are ~ and ~ it for [3] 
(version 3) and [I], respectively. For this reason, in the "best fit"-plot we do not fix 5 to 
its best fit value, but allow it to vary in its la-range. For the nu-plots 5 is allowed to 
vary in its ner-range. 

The last parameter which remains to be discussed here, is the smallest neutrino mass 
m . The strongest constraints on the absolute values of the neutrino masses comes from 
cosmology, where the sum of the light neutrino masses — in the form of the relic neutrino 
energy density Q u — is one of the parameters of the standard model of cosmology [8] . There 
is no unique consensus on the resulting upper bound on m v , however most constraints 
are of the order of < 0(1 eV) [8]. Therefore, we allow m to vary between zero 

and 0.3 eV. In the limit of m — > 0.3 eV, this implies an upper bound on the largest 
neutrino mass of 

I Atyi^ I 

maxm k <m + - — w 0.304 eV. (17) 

k 2m 

According to equation (Q, this directly translates to an upper bound of ~ 0.3 eV on the 
absolute values of the elements of the neutrino mass matrix. Since neutrinoless double beta 
decay has not yet been experimentally verified, there is no upper bound on the absolute 
neutrino mass scale from ((3f3)ov However, its non-observation gives upper bounds on 
|(-Mi/) ee |. The current bounds are of the order of 0, [TU] 

\{M v ) ee \ <0.4eV, (18) 
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la 2a 3a 

\{M v )ee\ (inv. spectrum) Forero et al. 1.52 x 1(T 2 1.36 x 10~ 2 1.14 x 1(T 2 

Fogli et al. 1.62 x 1(T 2 1.44 x 10" 2 1.24 x 10~ 2 
|(A4y)rr| (norm, spectrum) Forero et al. 

Fogli et al. 1.86 x 1(T 2 1.27 x 10" 2 



Table 1: Numerical results for the lower bounds on |(A^ i/ ) ee | (inverted spectrum) and 
| (.M,,).,-,- 1 (normal spectrum) in units of eV using the global fits of Forero et al. (version 
3) [3] and Fogli et al. [I]. The lower bounds for all other K-Mj,),^ are zero. 

which is comparable to the cosmological bound. Thus, unfortunately, the bound ( fl8|) will 
be fulfilled automatically in our numerical analysis, consequently giving no interesting 
physical insights. The bounds on mo coming from tritium decay are still of the order of 
rrio < 2 eV [8] and will therefore not constrain the results of our analysis. 

4 Results 

Lower bounds on K-Mi,)^: Using inequality (fTBl and the global fits of [31 H] shows 
that the only elements of K-M^)^! which have a non-trivial lower bound are | (^W^) ee | 
in the case of an inverted neutrino mass spectrum and |(-M;,) rr | in the case of normal 
ordering of the neutrino masses. The resulting bounds can be found in table [TJ 

How to read the plots: Before we discuss the resulting plots, we have to explain how 
they are to be read. This may best be done by means of an example. The properties we 
will outline in the following, hold for all plots in this paper. 

Figured] shows the allowed areas for |(A / ( I /) ee | vs. |(A^,/) Tr | in the case of a normal mass 
spectrum. The plot consists of three types of points, each one describing the boundary of 
an allowed area. The best fit area is bounded by black stars *, the la area is bounded by 
red triangles ▲, and finally the 3a area is bounded by blue points •. Note that also the 
plot axes comprise parts of the boundaries of the allowed areas, as can be seen in figure [TJ 
This feature is very pronounced for instance in figure [2J 

We present plots based on the results of the global fits of Forero et al. (versions 2 
and 3) [3] and Fogli et al. jl] in appendices lA j IB1 and IC j respectively. 

Plots of the smallest neutrino mass versus K-Mi/)^): Concerning the plots of 
K-Mj/)^! as a function of the smallest neutrino mass, an interesting point is whether they 
have changed since the original analysis of Merle and Rodejohann [5] in 2006. The main 
difference is that at the time when reference [5] was written, only an upper bound on 
the size of 6* 13 was known, which is the reason that [5] contains plots for different sizes of 
sin 2 26 , 13 . Comparing our plots to the ones of [5] with sin 2 26 l 1 3 = 0.1, we find that at the 
3a level the plots in [5] are still in good agreement also with the new results. However, 
for the best fit only the plots of the smallest neutrino mass versus |(-M^)ee| and K-Mj,)^! 
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Figure 1: Allowed ranges for |(A^ i/ ) ee | vs. KAl^rrl in the case of a normal mass spectrum, 
using the global fit of Forero et al. (version 3). The best fit area is bounded by black 
stars, the la area by red triangles, and the 3a area by blue points. 



are in good agreement with the plots obtained from the new data. The la regions are 
not indicated in [5]. 

Finally, from the plots we can readily read off the lower bounds on the K-Mj,)^ and 
compare them with the bounds found from evaluating inequality (fl3|) — see table [TJ We 
find full agreement, which provides a successful consistency check of our results. 

Plots of |(.M„) a/ g| versus \(M. u )a' In this section we will discuss some of the con- 
clusions one can draw from the correlation plots shown in appendices |A] [B] and O These 
plots can be used to test the viability of two texture zeros in M. v . The results obtained 
from the plots are in perfect agreement with the results of a recent numerical analysis 
provided in [IT] , and may — at 3a — be condensed to the fact that the seven two texture 
zeros originally presented in [12] are still viable with all of the three global fits [31 Sj. For 
further details we refer the reader to 

Let us continue by discussing those correlations which appear manifest at the 3a level. 
Going through all the plots, we find the correlations 
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Figure 2: Allowed ranges for |(-Mj,) er | vs. |(Aft,) Tr | in the case of a normal mass spectrum, 
using the global fit of Forero et al. (version 3). For more information cf. figure [T] and 
section @J 



most stringent. All these five correlations may be subsumed as follows: "If one matrix 
element is small, the other one must be large," as can, for example, be seen from figure HJ 
However, there are also matrix elements which appear to be totally uncorrelated. A good 
example for this case is |(A^ v )er| vs. |(A^ v ) Tr | in the case of a normal spectrum — see 
figure El 

While at the 3a level all plots produced using the three different global fits of Forero et 
al. and Fogli et al. agree very well, at the la level this situation changes completely. 
Instead of presenting a confusing list showing all differences, let us just mention the most 
important point which is the fact that the fit results of Fogli et al. no longer allow a 
vanishing matrix element |(.Mi,) rr | at la in the case of a normal neutrino mass spectrum, 
while the fit results of Forero et al. still do. This evidently also has strong consequences 
on all correlation plots including |(A4y) Tr |, a fine example being the plot of KA^)^! vs. 
K-M^tt-I, which is shown in figure |3j 

Let us finally discuss the sensitivity of the correlation plots to the data. As already 
pointed out, correlations which contain |(A^jy) rT | are particularly sensitive, due to the 
strong constraint on | {M. V ) TT \ in the case of a normal neutrino mass spectrum — see tabled! 
However, also for an inverted mass spectrum the plots based on the fits provided by 
Forero et al. (version 3) and Fogli et al. differ visibly. Other correlation plots which are 
quite sensitive to the data at la are: 
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Figure 3: Allowed ranges for |(.M„) MT | vs. |(.M„) rT | in the case of a normal mass spectrum, 
using the fit results of Fogli et al. For more information cf. figure [T] and section HI 
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A Plots based on Forero et al. (version 2) 

This appendix contains all 42 plots produced, based on version 2 of the global fit by Forero 
et al. [3]. 

The figures in the left and right columns correspond to normal and inverted mass 
ordering, respectively. The plots consist of three types of points, each one describing the 
boundary of an allowed area. The best fit area is bounded by black stars *, the la area 
is bounded by red triangles A, and finally the 3a area is bounded by blue points •. The 
same applies to the plots in appendices [B] and O 
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